Program package 
----------------------------------------------------------------------------------------------------
Program gic3d 
----------------------------------------------------------------------------------------------------
gic3d is a code to integrate GNSS velocity and InSAR LOS rate data for 3-D velocity field solution. 
To run the code, execute:
% gic3d < gic3d_b3.drv
----------------------------------------------------------------------------------------------------
Input files
----------------------------------------------------------------------------------------------------
gic3d_b3.drv                     ! driver file
v_swus.cmb                       ! GNSS station coordinate and horizontal velocity data file
v_swus.datv                      ! GNSS station coordinate and vertical velocity data file
subs.area1                       ! excluding area within which GNSS data are not used for InSAR ramp parameter estimation
crp_flt.dat13                    ! Coordinates of creep fault segments
d170.avg                         ! InSAR data input file #1 (data averaged in 0.02 x 0.02 degree grid)
a120.avg                         ! InSAR data input file #2 (data averaged in 0.02 x 0.02 degree grid)
a349.avg                         ! InSAR data input file #3 (data averaged in 0.02 x 0.02 degree grid)
d399.avg                         ! InSAR data input file #4 (data averaged in 0.02 x 0.02 degree grid)
d356.avg                         ! InSAR data input file #5 (data averaged in 0.02 x 0.02 degree grid)
d173.avg                         ! InSAR data input file #6 (data averaged in 0.02 x 0.02 degree grid)
d165.avg                         ! InSAR data input file #7 (data averaged in 0.02 x 0.02 degree grid)
d71.avg                          ! InSAR data input file #8 (data averaged in 0.02 x 0.02 degree grid)
a64.avg                          ! InSAR data input file #9 (data averaged in 0.02 x 0.02 degree grid)
a166.avg                         ! InSAR data input file #10 (data averaged in 0.02 x 0.02 degree grid)
d42.avg                          ! InSAR data input file #11 (data averaged in 0.02 x 0.02 degree grid)
d144.avg                         ! InSAR data input file #12 (data averaged in 0.02 x 0.02 degree grid)
d167.avg                         ! InSAR data input file #13 (data averaged in 0.02 x 0.02 degree grid)
d168.avg                         ! InSAR data input file #14 (data averaged in 0.02 x 0.02 degree grid)
d115.avg                         ! InSAR data input file #15 (data averaged in 0.02 x 0.02 degree grid)
a64wl.avg                        ! InSAR data input file #16 (data averaged in 0.02 x 0.02 degree grid)
d42wl.avg                        ! InSAR data input file #17 (data averaged in 0.02 x 0.02 degree grid)
d168f.avg                        ! InSAR data input file #18 (data averaged in 0.02 x 0.02 degree grid)
d169.avg                         ! InSAR data input file #19 (data averaged in 0.02 x 0.02 degree grid)
d167f.avg                        ! InSAR data input file #20 (data averaged in 0.02 x 0.02 degree grid)
-------------------------------------------------------------------------------------------------------
Output files
-------------------------------------------------------------------------------------------------------
gic3d_b3.out                      ! Integrated 3-D velocity solution
gic3d_b3.res                      ! InSAR data, model predicted, and Postfit residuals
gps_h_b3.out                      ! Interpolated horizontal velocities at GNSS sites 
gps_v_b3.out                      ! Interpolated vertical velocities at GNSS sites 
gps_b3.prd                        ! Interpolated 3-D regional GNSS velocities 
-------------------------------------------------------------------------------------------------------

-------------------------------------------------------------------------------------------------------
Program visr_n
-------------------------------------------------------------------------------------------------------
visr_n is a code to smooth a horizontal velocity field using an exponentially decay distance function. 
To run the code, execute: 
% visr_n < visr_n5a.drv
--------------------------------------------------------------------------------------------------------
Input files
--------------------------------------------------------------------------------------------------------
visr_n5b.drv                     ! driver file
gic3d_b3.out                     ! velocity data input
crp_flt.dat15                    ! Coordinates of creep fault segments
--------------------------------------------------------------------------------------------------------
Output file
--------------------------------------------------------------------------------------------------------
visr_n5b.out                     ! smoothed velocity data output
Entries:
Ve: east velocity (mm/yr)
Vn: north velocity (mm/yr)
See: east normal strain rate (micro-strain/yr)
Snn: north normal strain rate (micro-strain/yr)
Sen: shear strain rate (micro-radian/yr)
W: rotation rate (micro-radian/yr)
====================================================================================================
Data files
----------------------------------------------------------------------------------------------------
Data S1. GNSS velocity data used in this study. 
Data S2. Integrated velocity solution. 
Data S3. Strain rate field interpolated from velocity solution. 
====================================================================================================